n=[0.0 0.0;
   0.4 0.0;
   1.0 0.0;
   0.0 0.4;
   0.6 0.3;
   1.0 0.6;
   0.0 1.0;
   0.5 1.0;
   1.0 1.0];

u = zeros(size(n));
u(:,1) = 5*n(:,1) + n(:,2);
u(:,2) = -n(:,1);
sol = [n u]